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Abstract 



We demonstrate an alternative method for calculating the asymptotic 
behaviour of the discrete one-coin quantum walk on the infinite line, via 
the Jacobi polynomials that arise in the path integral representation. We 
calculate the asymtotics using a method that is significantly easier to use 
than the Darboux method. It also provides a single integral representation 
for the wavefunction that works over the full range of positions, n, includ- 
ing throughout the transitional range where the behaviour changes from 
oscillatory to exponential. Previous analyses of this system have run into 
difficulties in the transitional range, because the approximations on which 
they were based break down here. The fact that there are two different 
kinds of approach to this problem (Path Integral vs. Schrodinger wave 
mechanics) is ultimately a manifestation of the equivalence between the 
path-integral formulation of quantum mechanics and the original formu- 
lation developed in the 1920s. We also discuss how and why our approach 
is related to the two methods that have already been used to analyse these 
systems. 

1 Introduction 

The discrete quantum walk has been discussed in several recent papers 

The first authors to discuss the quantum random walk were Y. Aharonov, 
Davidovich and Zagury, in [2] where they described a very simple realization 
of the quantum random walk in quantum optics. Some further early results 
were due to Meyer, in ^J. He proved that for a discrete (unitary) quantum 
walk on the line to have non-trivial behaviour, its motion must be assisted by 
an additional "coin" degree of freedom which is conventionally taken to be two 
dimensional. This spin-like degree of freedom is sometimes called the chirality, 
and can take the values RIGHT and LEFT, or a superposition of these. Meyer 
therefore considered the wave function as a two component vector of amplitudes 
of the particle being at point n at time t. Let 



where we shall label the chirality of the top component LEFT and the bottom 
RIGHT. This paper is concerned with the dynamics of a test particle performing 
an unbiased quantum walk on the integer points on the line. At each time step 
the chirality of the particle evolves according to a unitary Hadamard transfor- 
mation 
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Figure 1: The discrete quantum walk on the line. The probability distribution is 
shown for a walk that started at the origin with its coin in the state |R) after it 
has evolved for 100 steps. The distribution is oscillatory between the two peaks 
and decays exponentially outside that range. The peaks move away from the 
origin with speed t/y/2, and the width of the peaks also decreases with time. By 
contrast, the classical random walk has a Gaussian distribution, which spreads 
with velocity ~ \fi. 



and then the particle moves according to its (new) chirality state. Therefore, 
the particle obeys the recursion relations 

* L (n, t + 1) = -^L(n + 1, t) + ^*R(« " 1= t) (4) 

* R (n, t + 1) = + 1, t) + ^*r(" - 1, *)■ (5) 

Meyer and subsequent authors have considered two approaches to the Hadamard 
walk, the Path Integral and Schrodinger approaches, which reflect two comple- 
mentary ways of formulating quantum mechanics |S] . We refer to the paper by 
Ambainis, Bach, Nayak, Vishwanath and Watrous |3] for proper definitions and 
references. We shall refine the asymptotic analysis of this paper. 

The behaviour of the Hadamard walk is very different from the classical ran- 
dom walk on the integer points on the real line. One way of understanding this 
is as a result of quantum interference. Destructive interference suppresses the 
probability amplitude in the vicinity of the origin, and constructive interference 
reinforces it away from the origin. The net effect of this is that the quantum 
walk spreads out much faster. Figure ^ shows the discrete quantum walk on the 
infinite line at t = 100. We have only plotted the distribution for even values of 
n in Figure^ If a walk's initial distribution has its support confined to a set of 
nodes which all have the same parity (all ns either even or odd) the support of 
the distribution will "tick" between different parities at each step. 
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The probability distribution for the quantum walk depends not only on the 
evolution law in J2J), but also on the initial conditions. Throughout this paper 
we will only consider walks that start at n = with the coin in the initial state 
|R). The convention that the walk starts in the position n = is inherited from 
the motivation for studying these systems as toy models of quantum algorithms: 
the computer is always started with its registers in the state |00 . . .). The choice 
of initial coin state was made because the Hadamard walk is an unbiased walk 
|12|. This means that even though some starting conditions result in an asym- 
metric probability distribution, we can always find some other starting state that 
will produce a walk with the opposite bias. The |R) starting state produces a 
distribution that is maximally biased to the right. Likewise, if we had chosen to 
start the walk in the state |L), this would have produced a distribution that was 
maximally biased to the left. This distribution is the exact mirror- image of that 
produced by starting in the state |R). Thus, reversing the starting condition 
just relabels n to —n. As the coin space is two-dimensional, we can now invoke 
the linearity of quantum mechanics, and note that we can obtain the behaviour 
for any initial condition |$) = a|R) + b|L) by forming the corresponding linear 
combination of the evolutions for the initial condition basis states |R), |L). 



2 Related Work 

We now very briefly describe the methods previous authors have used. They 
have so far followed two approaches to determine the limiting behaviour of the 
ip- functions as t — > oo. The translational invariancc of this problem means that 
it has a simple description in momentum space, and the Schrodinger approach 
relies on that fact. We will describe this in more detail in the section below. 
Beginning with the recursion relation (J1J, Nayak and Vishwanath showed 
that 

*(M) = (M fc )**(M) (6) 

where 

*(k,0) =^*(n,0)e ifc ™ (7) 

n 

and Mk is the matrix 

1 /_ p-ik P -ik\ 

M k = ^ ( J fe \ ik ) . (8) 

They diagonalize the matrix M k , finding the eigenvalues \\ = e luJk and 
Af = e^-^ where uj k = arcsin , € [— 7r/2, tt/2]. They then write 

the ^-functions in terms of the eigenvalues (A^) and (A|) and their associated 
eigenvectors and formally invert the original fourier transform to obtain the 
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closed form integral representations for the wavefunction, 
! + (-!)"+* r dk cosfc 



2 ,/_„. ^7T vi + cos 2 k 



1 + (-!)"+* r dk e lk 



feM ^'^" / ^ 71 = e-^W. (10) 

^ J- n 27r VI + cos k 

They then approximate these, using a combination of the method of stationary 
phase in one range, and integration by parts in the other. (Note that the Left- 
Right labelling convention in ^2] is the opposite to ours.) 

There is another approach based on the Path Integral formulation of quan- 
tum mechanics, that Ambainis et al. discuss. The V'-frinctions are expressed in 
terms of Jacobi polynomials (as in Lemma [3 below) of the form 

jW + ^'+«(0), j = (t n)/2 1 = (1 - 7 )t/2 - 1. (11) 

One may then derive the asymptotic behaviour of the ^-functions as t — > oo by 
determining the asymptotic behaviour of these Jacobi polynomials as m — > oo. 
This has been done in two ways. Ambainis et al. use the approach due to Chen 
and Ismail 7 , which employs the Darboux method. Here one begins with the 
Srivastava-Singhal generating function 



oo 

E 

3=0 



J^ +aj ' f3+bj \w)z j = (l + u)-T(l+v)- /3 [l + (l + a)u+(l + 6)v]- 1 (12) 



where u and v are defined to be power series in w that satisfy the equations 

u = -i(w + l)z(l + u)- A (l + v)-^ 1 (13) 

v= ™(u>- 1)^(1 + u)- A - 1 (l + v)~'\ (14) 

Chen and Ismail use the Darboux method to calculate the asymptotics of these 
Jacobi polynomials. This method starts from the idea that if f(z) = X) n >o a « z " 

then the radius of convergence is 1Z — linin^oo ( y^— r J when this limit exists. 
Suppose there is a comparison function g(z) — J2 n >o ^nZ n such that g(z) — f{z) 
has a larger radius of convergence than f(z), then b n — a n = 0(s~ n ) where s > 
7Z. If the asymptotic behaviour of b n is known, then since a n ~ b n then we know 
the asymptotic behaviour of a n . Chen and Ismail use the Srivastava-Singhal 
description of the generating function for Jacobi polynomials to determine its 
singularities on the radius of convergence and give comparison functions at each 
singularity to determine the asymptotic behaviour of the Jacobi polynomials. 

It is interesting to note that the reciprocals of these two singularities (when 
normalized by dividing by (V2) ) are the eigenvalues that arise in the Schrodinger 
method. 
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3 Overview and Results 



The other way to obtain the asymptotic behaviour of these Jacobi polynomials is 
to follow the method in Gawronkski and Shawyer's paper, |9]. This is the third 
way to analyse these systems, and will be the way that we follow in much of this 
paper. It is a refinement of the methods in the paper by Saff and Varga |14|. 
This uses the method of steepest-descents. We will outline this method below 
but for further details we recommend EH an d also the book by Olver J3j 
which describes both the steepest-descent method and the method of Darboux 
very clearly. 

The saddle-points that feature in this method are also related to the eigen- 
values that arise in the Schrodinger method. We will detail this relationship 
below. If C is a saddlepoint then exp(/i(£)) is an eigenvalue from the Schrodinger 
method, for a function h(Q that we define below. 

We now describe our results. We shall see that it is possible to obtain explicit 
asymptotic expansions that are uniformly convergent. The system displays two 
types of behaviour, with the transitions between the different behaviours gov- 
erned by a parameter, a = n/t. The behaviour changes qualitatively over three 
ranges, which are respectively < \a\ < 4^ — e, 4= — s < \a\ < + e, and 

■^j+e < \a\ < 1 — £, where e is a positive number. Our methods give error terms 

of the form 0(t~ N ~ 1 / 2 ) where N is some positive integer, and the expansion for 
each range holds uniformly. We stop at 0(t~ 5 / 2 ) since we have no application 
for the more precise estimates. 

The first range is arguably the most interesting. Here the asymptotic be- 
haviour of the ip- functions is oscillatory as t — ► oo. We will use a result by 
Gawronkski and Shawyer [U] who obtained the leading term and the factor 
1 + 0(i _1 ) for the error term for Jacobi polynomials. When this result is ap- 
plied to the ^-functions, we obtain a refinement of Theorem 2 of 3 who found 
the leading term. 

The second, boundary range is treated using the method of coalescing 
saddle-points as described in R. Wong's book |16| . A uniform asymptotic expan- 
sion is possible which involves the Airy function. This is also interesting because 
the ^-functions change from an oscillating, polynomially bounded asymptotic 
behaviour to an exponentially small behaviour as \a\ changes from below l/\/2 
to above l/v2- The calculation of the asymptotics for these polynomials is 
novel to the best of the authors' knowledge. The third range is the immediate 
vicinity of \a\ = 1/V2- As the behaviour of our integral representation is well 
understood, we can obtain uniformly convergent asymptotics thoughout this 
transitional range. 

The main interest of our results lies not in the more accurate asymptotic 
expansions for the ^-functions in the first and third ranges of a, although Am- 
bainis et al. did ask for a uniform method to do the asymptotics for these ranges. 
In fact, Gawronkski and Shawyer have already provided this method [5]. The 
fact that the uniform asymptotic behaviour for the Jacobi polynomials for the 
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transistion from polynomially bounded oscillating behaviour to exponential de- 
cay can be found in terms of the Airy function seems rather more interesting. 
It is somewhat surprising that this has not (to our knowledge) been published 
before for this family of polynomials: Wong's book 16 shows how to obtain the 
asymptotic behaviour in terms of Airy functions. Perhaps the reason for this 
omission is simply the previous lack of an application, which the quantum walk 
now provides. 



4 The Feynman Path Integral 

We begin with the Feynman path integral appproach following Meyer. We will 
represent the components of the vector-valued wavefunction as a normalized 
sum over signed paths. Meyer proved 

Lemma 1 (Meyer Let —n <t<n. Define I = The amplitudes of 

position n after t steps of the Hadamard walk are: 



Mn,t) = ^^[ t - l y-iy- (15) 
Mn,t) = ^Y,( l ~ 1 )( t ~ l )(-V t - S - 1 > (16) 



except for the endpoints where t = n, which have to be handled separately. 

(See for the endpoints. For a derivation of these formulae, see Appendix 
A of©.) 

We will follow the approach used by Ambainis et al. [3] and Meyer JT] but 
in greater detail to obtain the following lemma. The standard notation in for 
Jacobi polynomials is to write them as jb a '^\z), but to avoid confusion with 
a = n/t we will write these as jl u,v \z). (Note that in the following lemma, we 
have reintroduced the external phase that was omitted in [2] and 

Lemma 2 (Ambainis et al.[3 ). When n = t mod 2 and jl u ' v \z) denotes a 
Jacobi polynomial, 



^(n,t)(-i) (t - n)/2 

Also, 



[- ( |±a ) 2-™/2-i J ( t ^ n ) /2 _ i(0 ) ! when <n<t 

{-2 n / 2 - 1 jgl n f /2 _ 1 (0), when - t < n < 0. 

(17) 

' {-l) n+1 2- n / 2 - 1 J^+jl_ 1 (0), when 0<n<t 
( _ 1 y i +i 2 n/ 2 j(o + -«-i) (0)) whm _ f < n < q 

(18) 
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Remark: As Ambainis et al. [5j point out 

t — Ti 

M-n,t) = -Mn-2,t), M~n,t) = — — i> R (n,t). (19) 

t + n 

Proof of Lemma [21 We use two formulae from Abromowitz and Stegun 
to prove these results. The first is the Pfaff-Kummer transformation 3J (15.3.4 

of D) 

2 F 1 (a,b;c;z) = (1 - z)- a 2 F 1 (a,c ~ b;c; -^j). (20) 
The second is the representation of a Jacobi polynomial as a 2 F\, see 15.4.6 of 

m 

2 F x {-j, u + l + v+j;u+l;z) = f jj"^(l - 2z). (21) 

Now the first sum, say S, in Lemma^is 2 Fi(l — I, —(t — I); 1; —1) so by (121 )[l 

5 = 2( t -™)/ 2 - 1 3 Fi(l - (i - n)/2, -l+(t- n)/2 + n + 1 + 1; 1; 1/2). (22) 
Now by (EU 

we obtain 

s = 2 {t - n)/2 - 1 Jt L in-^- ( 23 ) 

This proves the first part of Lemma|21for n > 0. To derive the second part, for 
n < we see 

<S = a*i(-(* + n)/2, 1 - 1 + (n- t)/2; 1; -1) 
= 2 (*+")/ 2 2 ^(-(i + n)/2, i + („ +t )/ 2 -n-l; 1; 1/2) = 2^VVg^- 1 >(0). 

(24) 

The results for the first sum are now proved. 

The second sum, say T, can be treated in much the same way since 

T= (l-t) 2 .Fi(l-M-*+l;2;-l) = ^{l-t) a Fi(l-l,2-/+t-l;2;l/2) 
= 2 {t - n ^ 2 -\l - t) 2Fx(-{l - 1),2 + f - 1 + n; 2; 1/2) 

= (J-i)2(*-»)/ 2 - 1 J^ )/2 _ 1 (0)/I. (25) 

The result for nonnegative n now follows. If n is negative 

T= 2 F 1 (-(n + t)/2 + 1,1- l;2;-l)^- 
t + n 



2 

t + n 



2 (t+n)/2-i 2Fl (_(i + n )/ 2 _ 1, 2 + (n + t)/2 - 1-; 2; 1/2) 



2 



2 (t+«)/2- 1 aFl (_(( t + n)/2 - 1), 2 + ((n + t)/2 - 1) - t; 2; 1/2) 



t + n 

t + n" "(t+")/2 



2(*+™)/ 2 -ij(^) i(0) . (26 ) 
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This completes the proof of Lemma [21 

We consider first of all ipL(n, t), n > 0. (As was observed by Ambainis et 
al, there is a symmetry between positive and negative n.) We let 

m = (1 - a)t/2 - 1 (27) 

so that 

1 + a 2a ,„„. 

l + n = l + crf=- 1 to (28) 

1 — a 1 — a 



and 



4.1 The oscillatory range:|a| < l/y/2 — e 

For < a < 2~ x / 2 — e, < e < 2 -1 / 2 , we may use either the Chen- Ismail 
results or the Gawronkski-Shawyer [H] results. We use the results by Gawronkski 
and Shawyer since they have been proved to hold uniformly over this range of 
a. For the sake of consistency of notation and ease of reference we will state the 
Gawronkski and Shawyer result in a more restricted form than they obtained 
in their paper, but it is sufficient for our purposes. 

Gawronkski and Shawyer write, using the integral representation in equation 
(4.46) of Szego's book {T^, 

jiam+j.bm+p)^ = _L / e" 1 ^) g(()d( (30) 



r 



where 



ftK ^ 1,1 (^)) +< "°(^) +i,1 °(^) (31) 

and r is a contour circling the origin. These integrals are of the form 

I(m) = [ g(C)e- mh ^ d( (33) 



and as such can be approximated in the limit as to — > oo using the method of 
steepest descents [31 ED ^1- This relies the fact that in this limit (i.e., t — > oo) 
the only parts of the integrand that contribute significantly to the integral are 
those regions where the function in the exponent, h, has a maximum. This is 
because in the long time limit the exponential term in the integrand behaves 
more and more like <5-function(s) centred on point(s) where h is maximal. (Note 
that a stationary point of h is only a maximum along a given contour of integra- 
tion, r. This is because the stationary points of an analytic function can only 
be saddle-points, so whether they appear to be maxima or minima depends on 
the path taken through them.) 
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In order to make use of this phenomenon, we need to be able to assume that 
the imaginary part of h is approximately constant in the vicinity of these saddle- 
points, otherwise the integrand will oscillate unmanageably in the asymptotic 
limit. (We don't care if it oscillates wildly elsewhere, as the contribution from 
those regions will be negligibly small.) The way to achieve this is to choose the 
contour so that it passes through these saddle-points along the path of steepest 
descent for the real part of the exponent. 

We therefore choose the path T so it goes through the two saddle-points 
£ + ,C~ determined by 

h'(C) = 0( or (1 + a)( 2 - 2a( + 1 - a = Q). (34) 

Thus the imaginary part of h(Q is fixed which implies that the real part of h(C) 
has a maximum at the saddle-point. (The numbers e h ^ + \ e h ^ are the recip- 
rocals of the singularities found by Chen and Ismail using Darboux's method.) 
Gawronkski and Shawyer use a steepest descent contour which goes through 
the saddle-points at ( = ( + or £~ and the points £ = ±1. The contour must be 
modified slightly near the singularities at £ = ±1. The contours are leaf-shapes 
defined by 

3(MC)) = 3(MC + )). (35) 

So (for example) when a = 1/2 this is 

3?(ln(C 2 -l)-ln2-lnC+21n(l+C)) = 3(ln((C + ) 2 -l)-ln2-lnC + + 21n(l + C + )) 

(36) 

where 

+ _ a + Wl-2a 2 _ l + iV2 
1 + a 6 

Figure 12 is an example of the steepest descent curve for the oscillatory range of 
a. 

The Gawronkski and Shawyer result is stated in terms of the parameters 
A, B, £>, Ci, C 2 ,0, p{6) and £(0) defined as follows: 

2 

cos = - - a 9 SO ^r<0<TT, (38) 

1 — a 1 2 

A = Q, B = a, D = l-a 2 , (39) 

Gi=l, C 2 = \±^. (40) 
1 — a 

They then show that p{0) increases monotonically from to it as increases 
from to 7r or equivalently as \a\ increases from to I/a/2- Furthermore 

m = 9 + G ~ arctan { 1 + B 2B~ D tan0 / 2 )) > ( 41 ) 

./m t l + a(e (l + B 2 -D n .\\ 

m = - + - + —{-- arctan ( ^ tan 9/2 J J . (42) 

Here —tt/2 < arctan u < tt/2 for real u. 
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Figure 2: Steepest descent and ascent curves for the oscillatory range of a. The 
descent curve is the leaf-shape, and the two lines running off to ±00 are the 
steepest ascent curves. 



Theorem 1 (Gawronkski and Shawyer [9j). 

^(0+0.m,i±H + f^) ^ = + ^ ^ B ^ CQS 6/ / 2 )-l/2 x 

2am 1 + a 

(c^(l + B)((C 2 - 1)^4 + C 2 cos 2 0/2) 1/2 J 1_Q l ~" (sin(mp(0)+£(0))+O(l/m)) 

(43) 

as m (i.e., t) — » 00. 

Now, by direct substitution, 

/ 1 \l/2 

(7rm(l + -B) sin(0/2) cos(6»/2) _1/2 = 2 1 = (44) 



TriVl - 2cr 

Cz 1 (l + B)=l-a (45) 

((C a - l) 2 /4 + C 2 cos 2 (0/2)) 1/2 = -=± (46) 

V2(l - a) 

Thus from equations H3) , (011 , (H3) and JH3), 

1/2 



j(0,2q/(1-q) + (1+ Q )/(1-q)) ( ' n = 2 f 1 "\ 

V ' \ttWI - 2a 2 / 



x 



2 (2a((l-a)t-l))/(l-a) + (l+a)/(l-a))/2 g j n _ _ + ^ (l +0 (l/t)). 

(47) 
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We have that 

a 1 + a at 1 , . 



so 



Ji' — ^CO) = * _ 2 •Bin((l-a)t/2-l)p(g)+^(g))(l + Q(l/t)). 

(49) 

Now it follows that for — t < n < t, 
Mn,t) = —jJ2=-j- 2 sin((l - a)t/2 - l)p{9) + + 0(1/*)) (50) 



where 



... S 2a fa, 

p(9) = arctan , (51) 

FW I- a I- a \VT^ ' 

a(g) = -7— — + tt/4 - i±£ arctan f ) . (52) 



1 — a 1 — a \ 



a 



When we consider ipft(n,t), < n < t we find that the power of 2 does 
not change. (Gawronkski and Shawyer use the symbol a as a dummy variable 
in their Jacobi polynomials. To avoid confusion, we will call this parameter 
k.) For ip^(n,t), k was set to be zero. To do the calculation for i/jn(n,t) we 
need to set K = 1. Note also that we must reset r (which they call (3) to be 
t = b = 2a/ (1 — a). Thus p{9) does not change. However, £(0) does change as 
we will specify below. We (as Gawronkski and Shawyer do) use 

km ( 1 + A2 -f + D t a n e -) (53) 
2 A^o+ I 2A 2 / V ' 



to obtain 



2 VVT-Q' 



Cr(^) = Cl(^) - J + arctan ( -=== ) . (54) 



Thus we have: 



Theorem 2. Let e > be any constant and a be in the interval (-1/V2 + 
e, l/v2 — e). Then as t — > 00 we have uniformly for < \a\ < l/v2 — e, 

Mn, t) = Q±fj) 1/2 f sin ((^ - DpW + few) (1 + 0(1 A)) 



(55) 



Mn, t) = —^UL=— sin ({^t - l)p(g) + fe (g)) (1 + 0(l/t)) (56) 

with p, £l and fe defined by equations 1)5 1[), (|52() and (I54|l . iJere 

6» = arccos(-a 2 /(l-Q ! 2 )), ^ < < tt. (57) 



2 
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Remark: The term of the form l/yl — 2a 2 gives us forewarning that this 
term is going to become very large when \a\ = l/y2. This is consistent with the 
graph in figure^] In fact this term actually diverges at this value of a, but this 
is a symptom of the breakdown of this approximation, which is why we included 
the £ in the statement of the theorem. In this transitional range, Theorem [3] 
below is the appropriate form to use. 

Our Theorem |21 agrees with Theorem 2 of Ambainis et ai, as expected, and 
it also gives an estimate for the error term. We skip the proof that the answers 
are identical for the -^-functions as the probability function p{n 1 t) defined by 

P (n,t)=iPl(n,t)+iJjl(n,t) (58) 

is more interesting and we will show that for this function and its moments, our 
answers are identical. We will use the identity sin 2 A = (1 — cos2^4)/2. Recall 
that a — n/t, where we temporarily think of t as fixed and let n vary. We may 
then use B(n,t) to denote a bounded function with bounded derivatives and 
A(n, t) for a function such that it and its derivatives are bounded away from 0, 
thus for ipYi(n,t) we have 

^"'(Sl"' ^ (60) 

and for ijjh(n,t) we have instead 

B(n,t) = (62) 

This enables us to write: 

B(n, t) cos(A(n, t)t) da = 



r 2 , x d ( sin(A(n,t)t)\ , , , 

L B(n>t) ^{ tA'( n ,t) ) d " + °(V*) = 
B{n,t)s\n{A{n,t)t)Y 2M p M B'(n, t) sin(A(n, t)t) 



A'(n,t)t 



n 1 ,t 1 



„„„ A.(n, t)t «»+0(V*) = 0( 1/t ) 

(63) 



by a simple integration by parts, where a\ = n\/t\, 012 — "2/^2- If we write 
p(a) = tp(n, t) (following Ambainis et al.) and note that the O-term is uniform 
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for — 1/V2 + e < a < \/\/2 — e, then that observation and our Theorem [21 give 
us that 



ft 2 



provided \a\ < l/v2 — e. Note for the quantum walk p(a) is when n and t 
have unequal parity so for the quantum walk we have 



ft 2 



1 



ftl 



p(a) da=- - 7-7== da + 0{l/t) (65) 

7T J ai (l-a)Vl-2a 2 



To confirm that we have a probability distribution, we must verify that the 
function integrates to 1 : 

-e+i/V2 l 1—E+1/V2 l 

p{a) da ~ — / -== da (66) 

-1/V2 ' f Js-i/V2 (1 - a)Vl - 2a 2 

i /■- E+1/V2 i 

da. (67) 



W e -i/V2 (l-a 2 )Vl-2a 2 



If we let cosfc = —a/y/l — a 2 , and sinfc = vl — 2a 2 /\/l — a 2 , then we can 
write 

— — = ; 68 

da (l-a 2 )Vl-2a 2 

so 

1 /-e+l/\/2 1 ,«(-e) 

- / p(a) da = - dk=l- 2tt- 1 5(s) (69) 

K Je-1/V2 n J-n+6(e) 

where lim e ^n S(e) = 0, as required. The correction term appears because we 
have only performed the integration over the oscillatory range of the probability 
function, as this is where it has almost all of its support. 

We can now write down the m-th moment of the distribution: 



1 



1/V2 

da. (70) 



7T J-1/V2 (1 - a)Vl - 2a 2 
Thus the first moment is 

1 f 1 ^ a 2 -l , 1 f 1/V ^ 1 

da H — / , da 



W-1/V5 (l-a 2 )Vl-2a 2 tt J-1/V2 (1 - a)Vl - 2a 2 

= -^/ /1 o 2 da + l = l-VV2. (71) 

^ 7-i/^2 Vl - 2a 2 

The second moment can be seen to be also equal to 1-1/^/2. 
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4.2 The exponential range: |a| > 2/y/2 + e 

We now consider the range a > 1/V2 + e where < e < 1 - l/\/2. The 
Gawronkski and Shawyer results are an extension and refinement of the results 
of Saff and Varga ^3]- We state the results as Saff and Varga do. They write 

jjhr,3om/(l-a)+#( ) = y e mh ^g(t)d( (72) 

where 

2a 

h(() = ln(C 2 - 1) - In 2 - InC + ln(l + C), (73) 

1 — a 

fl(C) = 2^(l + /, (l-0 7 J- (74) 
They choose ( to be the saddle-point £~ , where 



1 + a 

Using the saddlepoint method they derive 
Theorem 3 (Saff and Varga [Tip. 



r = q-V2^-l (?5) 



J%<*™/V+°)+f>(0) ~ zexp[m/ l (r)] 5 (D (^lj^y) 1/2 (76) 



where 



Now 



x a - V2a 2 - 1 a + V2a 2 - 1 a 2 + 2\/2a 2 - 1 

C~C - 7- 1 ; 2 ' ( 78 ) 

1 + a 1 — a 1 — a 

l + C^^ 2 "'^^. (79) 

1 + a K ' 

Then according to Saff and Varga, 

/ \ m / \ 2am/(l-a) 

j( 7 ,2am/(i-a)+/3( \ _ [ a + V 2a 2 - 1 \ / 1 + 2a - V2a 2 - 1 \ 



1 — a 2 I \ 1 + a 

/3 / , * \ 1 



1 / 1 + 2a - V2a 2 - 1 \ / 1 + V2a 2 - 1 \ 1 



2tt\ 1 + q \ 1 + a / a - \/2a 2 - 1 V m/l "(0 



1/2 



(80) 

where h (£) is defined in (|7T|> . 
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Gawronkski and Shawyer show that the ~ symbol can be replaced by 1 + 
0(l/t) and that this expansion holds uniformly for a <G [l/v2 + e, 1 — el- 
Note that when a — » 1/ -\/2, 

J (7,2 Om /(i- Q )+(J)( ) 1/m ^ (ll^L} 1 " 2 =2 Q */ 2 = 2"/ 2 (81) 



so the asymptotic estimate above metamorphoses into the asymptotic estimate 
in Theorem [3 We shall state our results for positive a as those for negative a 
follow immediately. These results refine the estimates of Ambainis et al. [3]. 

Theorem 4. // a > 1/ y2 + e, then uniformly for a G [1/ \/2 + e, 1 — e] , 

Mn,t) = C R ^B t (a)(l + 0(l/t)) (82) 

Mn,t) = C L -^B t (a)(l + 0(l/t)) (83) 
where 5(a), Cl, Cr are defined in the asymptotic expansion of 

2 ~" /2 ~ 1 47' 2Qm/(1 " a)+/3) (0) (84) 

following from the above. (Note that for Cl, 7 = and (3 = (1 + a)/(l — a) 
ly/iiZe /or Cr, 7 = 1, /3 = 2a/(l — a) m i/ie above.) Thus 



(85) 




(86) 



(1 - a) (2a + 1 - \/2 



2 



a(V2a 2 -!-(! + 2a))(l + af \ 
(l + a(l-a + V2a 2 -l)) 2 ~) 



(87) 



to reconstruct the complete wavefunction. The asymptotics for a in the range 
— 1 + e < a < — 1/V2 — £ follow from the spatial symmetry between —n and n. 

Figure |3] shows the steepest descent curves for this range. 

Remark: Saff and Varga show that B(a) is a decreasing function of a for 
a 6 [1/V2 + £, 1 — e]. Thus the functions decrease exponentially in t in this 
range. 
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Figure 3: Steepest descent and ascent curves for \a\ > l/y/2. The teardrop 
shape and the line from £ + are both steepest descent curves, but only the 
contour that surrounds the origin gives a valid integral. The remainder of the 
lines are steepest ascent curves. 



4.3 The transitional range: 2 x l 2 — e < a < 2 x l 2 + e 

We now consider the range l/y/2-e < a < l/y/2 + e. Theorems and H exhibit 
an oscillating sine term times i -1 / 2 for \a\ < l/v2 — e and an exponentially 
small estimate for \a\ > 1/V% + £• Saff and Varga ^1] show that if a = \j\pl 
then 

J^~(0) ~ ^ (- p ^_ TJ ) ,/3 «-/ ar(1/3) 

(See also Ambainis-et al, ,3 .) The asymptotic behaviour is therefore qualita- 
tively different in the three cases |or| < l/y/2 — e, a = 1/V2, and \a\ > 
l/\/2 + e. 

We now apply the analysis for coalescing saddle-points as described in R. 
Wong's book |16|. to get a uniform asymptotic expansion for a range of a extend- 
ing a positive distance on each side of a — 1/V2. The asymptotic behaviour 
is described in terms of the Airy function, Ai(x). These functions cannot be 
written out explicitly, but (like Bessel functions) they can be written in terms 
of integral representations that are well understood. It suffices for us to know 
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that as x ~ * — oo the behaviour of these functions is oscillatory, 

Ai(-x) ~ -L cos(|x 3/2 - tt/A)/x 1/ \ (89) 

V7T 3 

Ai'(-x) ~ ^ sin(^x 3/2 - tt/4) (90) 

but that when a; — > oo they behave like 

Ai{x) ~ ^™ 1/4 exp (~^ 3/2 ) , (91) 

^'(^--^^exp^-Hx 3 / 2 ) (92) 

which decreases faster than any power of x. 

Recall that we have the two saddle-points £~, ( + = (a=p v2c? — 1)/(1 + a). 
Following Wong we define the variables C and rj by 

C 3/2 = ~(/>(C + )-MC)) (93) 



or 2C 1/2 = C + - C" = 2iy/l - 2a 2 /(l + a), (94) 
V=l(h(0 + h((+)) (95) 

We must choose (~ , £ + as we have because we want the range of a > 1/V2 
to give positive C so that there is exponential decay as C ~ * 00 • Note also (as 
Saff and Varga point out) that if a = l/y2 then 

h'(a) = h"(a) = (96) 

/T(a) = * + " (97) 

v ; 24a 2 (l-a) 3 v ; 

Thus 

h(l/V2 + 5) = - l + a 5 3 + 0(**) (98) 
24a z (l — a) 

so 

mh(a) = o(l) if 6 = o(m~ 1/3 ). (99) 

Also /i(C + ) ~ MC~) so r ? ~ ^(C~) f° r this range of a. 

Suppose for now that equation 4.31 in chapter VII of Wong ^|)] holds, which 
in our notation becomes 

J^Wd-H^o) = 2nie™> ^ff oo + ^P&o) (l + 0(m- 1/3 )) 

(100) 
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where ao and 60 are independent of m. (In Wong's book, our m is his A and our 
£ is his t. Our Jacobi polynomial is his /(A; a), but his a is something else.) To 
obtain an expression for -0l we need to set 7 = and 

/3=i±^ (101) 
1 — a 

as in Theorem El Likewise, for t/jr, we will need to set 7=1 and 

/3 = T^-- (102) 
1 — a 

The previous argument shows that the term mr\ — o(l) if <5 = o(m -1 / 3 ) and so 
the e mri term is asymptotic to 1. Note furthermore that if a = 1/V2 + 6 then 
the definition of £ implies that 

C = -±PK = 6 + Ojfi). (103) 

1 + a 2 (1 + V2) 2 



Thus 

m 



2/3 C = l^^^/sj + o^). (104) 

(l + \/2) 2 V ^ V j 

When we use this estimate in the asymptotic behaviour of Ai(m 2 ^ 3 Q and 
we see that the Airy functions give terms superpolynomially small if 
5 > m~ 2 / 3+c . If S — mT' 1 , 77 < 2/3 however we see that the ^-functions are only 
polynomially small. A similar argument shows that if \a\ < 2~ 1//2 — m~ 2 / 3+e 
then the behaviour of the tp is oscillatory. 

We can therefore write that the transition from oscillatory behaviour bounded 
below by a power of m or t to bounded above by a superpolynomially small func- 
tion occurs when a varies by 0{t- 2 / 3 ) from l/y/2. 

We may use equation 4.31 in chapter VII of Wong 16] if the transformation 

h(Q = u 3 /3 - (u + r] (105) 

is single-valued on the contour of integration. We may use the Gawronkski- 
Shawyer contour of integration (which Saff and Varga also use). The transfor- 
mation will be one-to-one if on the path of integration, 

^ ^ and % ± 0. (106) 

Now 

i = im (107) 

and the bottom derivative is only zero at the saddle-points. The only saddle- 
points are at (~ and £ + ■ Wong shows that the choice of C~ , C + we used implies 
that 4^ ^ 0. The only place the numerator is zero is at £~ , C + an d the denomi- 
nator is analytic and therefore bounded on the contour of integration. We find 
from Wong's book that 
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Theorem 5. There is a positive e so that uniformly for some clq and bo (defined 
below) 

j(7,W(i-<*)+/3) (0 ) = 2me~""i ( ^ffi^ o + ^|ff^ o) (1+0(™~ 1/3 )). 

(108) 

When n = t mod 2 and jl w ' v ^ (z) denotes a Jacobi polynomial, 

ipL(n,t) = 2-"/ 2 - 1 ji°'^' + ^ f (0), w/ie«0<n<t, (109) 

where we can use the Remark following lemma\^to obtain the wavefunction for 
negative n. Also, 

^H(n,t) = (l^)2- n / 2 - 1 J ( n 1 l '^ + ^\0), when0<n<t, (110) 
1 — a 

and use the symmetry properties to obtain the other half as before. 

So this integral representation is valid for all values of— l + e<a< 1 — e. 
In this theorem, 

Mu) = g(0§ 1 (in) 

and 

ao = llMC 1/2 ) + M-C 1/2 )}, (112) 

bo = llMC 1/2 )-M-C 1/2 )}- (113) 

Remark: Theorem [S] gives the asymptotic behaviour in a rather convoluted 
way. It is really only useful for \a\ very near to l/v2, where the behaviour 
undergoes a qualitative change from oscillatory to exponential decay. For other 
values of a, Theorems and 01 give much simpler expressions for the wavefunc- 
tion. 

Figure 01 shows the steepest descent curve when the two saddle-points coa- 
lesce. 

4.4 Error bounds for the method of steepest descents. 

We now discuss very briefly the error-terms in our results so far. We restrict 
our attention to ipi,{n,t) and Theorem |21 above, but similar comments apply to 
the other results in this paper and tpn(n,t). In order to apply Theorem 7.1 of 
Olver we expand h(() and g(£) in powers of £ — £ + : 

g(0 = 5(c + )+E c ^-c + r (us) 
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1 



Figure 4: Steepest descent and ascent curves when the saddle-points coalesce. 
The tear-drop shape and the line from comprise the steepest descent curve, 
and the two lines out to ±00, together with the line to the origin are the steepest 
ascent curve. 



and write 

/ 00 \ 00 

e £r= 3 Mc-c+r U(c+) + ^ Cr (c-c + r =^a s (C-C + ) s - (116) 



s=0 



The steepest descent contour that Gawronkski and Shawyer use naturally 
separates into two pieces, a piece Ti above the real axis and a piece T2 below 
the real axis. Notice that Ti is the mirror image of T2 in the real axis. With 
this notation, Theorem 7.1 in Olver's book says 

J r e-™^ 5(C) ^ 2e -™MC + )g r(s + 1/2) _|Ei_. (117) 

Here 

«o = jdyr* (U8) 

r „ 2b! 11 fb{h m ) 2 h""\ \ 1 ,„,„, 

when <?, /i and their derivatives are evaluated at Q — C + - The asymptotic expan- 
sion of the integral over T2 is the complex conjugate of the integral over Ti. It 
is therefore possible to derive complete asymptotic expansions for tpi,(n,t) and 
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tpn(n,t) and certainly Gawronkski and Shawyer were aware of this. The s = 
term is the asymptotic formula of Theorem [21 above. Since we do not need more 
precision than in Theorem[21for this application, we do not pursue this further. 

Olver then explains how to derive numerical estimates for the error term. 
This is somewhat complicated so we refer the reader to section 10 of chapter 4 
of Olver's book. One needs to compute the maximum of certain quantities on 
the contour of integration. Since we do not have applications of such bounds 
we will not pursue that here. Olver also shows how to derive explicit numerical 
error bounds when the asymptotic expansion immediately above is truncated 
at any value of s. 



5 The Schrodinger Approach 

Nayak and Vishwanath ^5] start from the recursion relations in (J2J and use the 
Fourier transform 

V(k,t) = ^^{n,ty kn (120) 

n 

where ij;(n,t) is defined by and obtain 

*(fc, t + 1) = M k ^{k, t) (121) 

when 

V2M k =^ k 6 ~f). (122) 

The eigenvalues of M k are 

X 1 = e iUk , A 2 = e i(7r-w fc ) ( 123 ) 

where ujk £ [— 7r/2, 7r/2] and satisfies sino;& = sin 
It follows that 

* k (t) = M${k,0). (124) 
They deduce from this that 

~ . . If cosfc \ ■ , , (—1)* / cosfc \ ■ , t 

* L (M) = - 1 + ; e tuJkt + - — - 1 - - e~ Wkt 

2 V vi + cos 2 fcy 2 V vi + cos 2 fcy 

(125) 



9 R (k, t) = (e*"» - (-Ife-^) . (126) 

Z\L ~\~ cos /c 

Formally inverting the original Fourier transform (using Cauchy's integral for- 
mula) and some ingenious manipulations produce 

j- , 2 = e-^-^ (127) 

-7T 2n V 1 + cos 2 k 

*a(n, t) = f £ ( 1 + , C ° Sfc 2; ) e^*-*») (128) 

J-tt 27r V V 1 + COS fc/ 
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where w fc = arcsin y^ff) € [— 7r/2, 7r/2]. 

They then apply the method of stationary phase to obtain a weaker version of 
Theorem|21above, and integration by parts to show that the wavefunction decays 
superpoly normally fast for |a| > l/y/2, which gives them a much weaker version 
of Theorem0](Ambainis et al. show that this decay is exponential, but they were 
unable to obtain uniform asymptotics). Of course both approaches consider the 
same functions t) and 'J'r^, t), the differences are just the representation 

of the generating functions and the choice of contour of integration, as we will 
discuss below. 

If \a\ < -j- each eigenvalue is minus the complex conjugate of the other, so 
the ^-functions have an oscillating factor. When we find the stationary points 
of the phase, we obtain an equation for k at the critical points, k a . This is 

cosk a = -f^=. (129) 
yl-a 

When \a\ < l/y/2 this has solutions which are real and distinct. The two 
solutions merge at a = l/y/2, and then become complex. When a is outside the 
range \a\ < l/y/2 the phase has no stationary point on the real axis. We have 
been unable to find a method for approximating these integrals. It is worth 
noting that (|127|l are themselves integral representations of Jacobi polynomials 
as a function of its parameters. 

The exponentially decaying solutions are counter-intuitive in other ways. 
As we mentioned above, for this case k a is complex, so instead of seeing the 
oscillatory behaviour we might be expecting, instead the wavefunction decays 
within an exponential envelope. This is rather like the phenomenon of evanes- 
cent waves. These can also occur classically: consider an electromagnetic wave 
incident on the surface of a conductor. These waves cannot propagate in con- 
ductors, as the latter will not sustain charge gradients. However, the wave does 
impinge a finite distance into the conductor (the "skin depth") over which its 
amplitude decays exponentially. Mathematically, this is equivalent to a com- 
plex wave-number. Evanescence can also occur in quantum mechanics, typically 
in regions that are classically forbidden to the particle. Strangely, both these 
scenarios involve the presence of some kind of barrier, but no such barrier is 
present for the quantum walk. These regions are not classically forbidden to 
the particle, it's just very unlikely to be there. 

5.1 The relationship between the two approaches 

In the paper by McClure and Wong [TJJ the authors show that the methods of 
stationary phase can be reduced to the method of steepest descent under quite 
general conditions, as the same results can be obtained from either method, 
with exactly the same convergence properties. We can see this intuitively as 
follows. Using steepest descents, we have the integral 

-1- / g(()e- mh ^ d(. (130) 
2ni J r 
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If the contour T goes through the saddlepoint ( + we can deform T to the curve 
ICI = IC + I = r or C = re 10 , for some dummy variable 9. This produces an integral 

— / e- m}{reZe) g{re ie )re t(> dd. (131) 

27r J -IT 

If we separate / into its real and imaginary parts, then we obtain 

_L J* e -rnmf{re iS ))+i*U{re*°))) g ( re i6y e ie ^ ^ 

For a quantum system undergoing unitary evolution, 3?(/) = and since r 
is a constant, this integral is of the form required for the stationary phase 
approximation. So we will write /j,(8) = Q(/(re* e )). 

The stationary points are those for which //(#)' = but since 

we see that the stationary points are identical to the saddle-points. With the 
method of steepest descents the integrand has a very small absolute value away 
from the saddle-point. By contrast, in the method of stationary phase the 
oscillations of the kernel /i become arbitrarily rapid away from the stationary 
point, and so self-cancel so long as g{re l9 ) is sufficiently smooth. (Readers 
requiring further details are referred to the lucid exposition in 



6 Conclusion 

We have developed a new way of analysing the discrete quantum walk on the 
infinite line in terms of Airy functions, which has the advantage of being able 
to handle the dramatic changes in the asymptotic behaviour of this system in a 
uniform manner. We have also probed the mathematical relationship between 
the path-integral and Schrodinger approaches to solving this problem. Previous 
authors have found the methods of integration by parts and stationary phase 
to be problematic over some parts of the range of a. By contrast, the method 
of steepest descents yields a unified treatment of the system. 
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